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A discrete chemotactic predator-prey model is proposed in which the prey secrets a diffusing 
chemical which is sensed by the predator and vice versa. Two dynamical states corresponding to 
catching and escaping are identified and it is shown that steady hunting is unstable. For the escape 
process, the predator-prey distance is diffusive for short times but exhibits a transient subdiffusive 
behavior which scales as a power law t 1 ^ 3 with time t and ultimately crosses over to diffusion again. 
This allows to classify the motility and dynamics of various predatory bacteria and phagocytes. 
In particular, there is a distinct region in the parameter space where they prove to be infallible 
predators. 
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Phagocytes or predatory microbes are hunting their 
prey by chemotaxis [T], Q , i.e. they sense the concentra- 
tion of a chemical which is secreted by the prey and is 
diffusing through the solution Q . Typically the predator 
moves along the steepest gradient of the chemical concen- 
tration to ultimately find its ejection source. Likewise the 
prey (for example another microbe) "smells" a secreted 
chemical from the advancing predator and tries to escape 
by moving along in the opposite direction of its max- 
imal gradient. This chemotactically coupled predator- 
prey systems are relevant for many biological microor- 
ganisms. In fact, there are many examples of biologi- 
cal relevance for chemotactically coupled predator-prey 
microorganisms. To name just a few, common micro- 
bial predators and phagocytes are Bdellovibrio [H-dj], P- 
aeruginosa Q, D. discoideum [1, lymphocytes [l(| 
and M. xanthus [Til fl^] ■ 

Previous theoretical investigations have focussed 
on spatiotemporal pattern formation in predator-prey 
colonies [HI, [Ti| which are typically described by nonlin- 
ear reaction-diffusion equations [TEI, [lj| . While the lat- 
ter approaches involve a coarse-grained continuum mod- 
elling, there are much less model studies on individual 
microorganisms. A discrete swarming model of individ- 
ual self-propelled particles for bacterial colonies has been 



proposed by Csirok ct al 13 1. This was elaborated re- 
cently by Romanczuk et al 17] based on a related individ- 
ual model of Schweitzer and Schimansky-Geier [l8j]. Fi- 
nally individual autochemotactic models have been stud- 
ied where the microbe follows its own diffusing secretion 
[l9l - l22l |. In all of these individual models there is no 
predator involved, apart from a recent study [23j which 
addressed a lattice model with no chemicals involved. 

Here we propose a discrete model which describes both 
the predator and the prey individually and contains ex- 
plicitly the diffusion of the two chemicals secreted by the 
predator and the prey together with the Brownian mo- 
tion of the latter. The deterministic (fluctuation-free) 
model is analyzed analytically and by numerical solu- 
tion which is supplemented by Brownian dynamics com- 
puter simulations at finite temperature. Depending on 



the model parameters and the initial distance between 
predator and prey, two different dynamical processes are 
identified which correspond to catching and escaping, and 
an unstable steady hunting. By analytical treatment, var- 
ious scaling laws are extracted characterizing and delin- 
eating the different regimes. In the absence of noise, the 
mean-square distance between predator and prey scales 
with different exponents a as a function of time with a 
subdiffusive anomalous exponent a — 2/3 for escaping 
24] and a ballistic behavior a — for steady hunting. 
Brownian motion leads to ultimate diffusion (a — 1) such 
that these exponents are transient. Catching is accom- 
panied by a scaling form of \t — t cap \ a , with a — 2/3, t cap 
being the capture time. In principle, our results allow 
to map and classify different biological systems into the 
different regimes. 

In our discrete predator-prey model, the predator is 
at position ri(t) at time t, hunting the prey which is at 
position r 2 (i) and trying to escape. The concentration 
field of the chemical secreted by the predator (prey) at 
a constant ejection rate Ai( 2 ) is denoted by Ci( 2 )(r,i). 
Assuming a gradient sensing scenario for each microbe 
in response to the chemical secreted by the other, and 
taking into account effective stochastic fluctuations that 
are associated with the non-equilibrium self-propulsion 
mechanism of each, the overdamped equations of motion 
for the predator and the prey, respectively reads 



71*1 = +KiVc 2 (ri,t) + -q^t) 
72*2 = -K 2 Vci(r 2 ,i) + f] 2 (t) 



(1) 

(2) 



Here, 7i( 2 ), ^1(2) an d *7i( 2 )(£) are the damping constant 
in the medium, the chemical coupling constant (a mea- 
sure of the gradient sensing strength) and the effective 
noise vector associated with stochastic self-propulsion of 
the predator (prey), respectively. The first term on the 
right, in both equations, models the systematic contri- 
bution of chemotactic response through simple gradient 
sensing. We take both ki >2 > 0, so that c 2 acts as a 
chemoattractant for the predator, while ci is a chemore- 
pellant for the prey. We model rj i (t) as a Gaussian white 
noise: (r] t (t)) = 0, (Vi»{t)Vj»{t')) = 2 ll ^ 1 8 l] 5^8{t~t l ). 
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The Greek indices refer to spatial components, while the 
Roman indices are reserved for the microbe's attributes 
(i = 1: predator; 2: prey). Here, /3 corresponds to an in- 
verse effective temperature associated with the stochas- 
tic fluctuations, such that Di = 1/(^/3) is the non- 
chemotactic diffusion constant of the microbe concerned. 
Hydrodynamic interaction between the microbes is ne- 
glected. 

The diffusion equation of each of the chemicals reads 



da{v,t) 

at 



+u i (r,t)-Vc i (r J t) = D cl V 2 c l (r,t)+X i 5[r-r l {t)} 

(3) 

where D C i is the diffusivity of the corresponding chemi- 
cal, and we have assumed each microbe as a point-source 
emitter. Uj(r, i) is the advective flow-field set in the 
medium due to the motion of the microbe. In practice, 
a microbe has a typical mesoscopic size 'a'. Expressing 
lengths in units of a: r — > r" = r/a, and time in units of 
t = a 2 /D ci : t—tt" = t/ro, Eq.(3) reduces to 



d_ 

dt" 



auj(r,t) 

D ci 



7" 2 



a 2 \i 



,(r,t) = -fij[r-r i (*)] 
(4) 

Simple spatial gradient sensing microbes are known to 
move slowly with velocity on the order v ~ 10 ~ 2 — 
10 /zm/s, and are typically of size a ~ I — 10/xm. Chem- 
icals are secreted typically at Aj ~ I0 3 molecules/s and 
diffuses at D ci ~ I0 2 — I0 3 /im 2 /s. Under such practical 
situations, noting that the magnitude of the flow-field 
| can be at most on the order of v, we have the dimen- 
sionless factor av/D C i ~ 10~ 5 — 10~ 2 . This makes the 
advective term negligible, and reverting to the original 
variables: r" — > r and t" — > t, Eq.(4) simplifies to 



dcj(r,t) 

at 



D cz V 2 Cl (r,t) = Xi6[r-ri(t)] 



(5) 



Fast moving microorganisms with v ~ 10 — 10 2 /ito/s (ca- 
pable of producing appreciable advection in the medium) 
are known to chemotax by 'temporal sensing' mechanism 
of the chemical gradient [l^ which we do not address 
here. Also, the detailed influence of the microbe's shape 
and distribution of chemical sensors on the spatial gradi- 
ent sensing itself are not considered in our model. 

For an unconfined space in three dimensions, the 
Green's function solution to Eq.(5) yields 



Cj(r,i) = A, 



dt' - 



exp 



4D ct \t-t'\ 



Brownian dynamics of the predator-prey system is imple- 
mented to simulate the chemotactic motion (eqs.(l) and 
(2)), and using Eq.(6) to calculate the spatial gradient of 
the chemical concentration. We measured time in units 
of A^ 1 , lengths in units of Iq = O.ly/D^/X^ an d energy 
in units of eo = 72^2- Real estimate yields K% ~ 10 3 eo'o 
for Dictyostelium [9( moving at 0.2/im/s up cAMP gra- 
dient of 0.01nM//im, secreted at 10 3 molecules/s with 




FIG. 1: (color online). (Left) Chemotactic Chase: A predator 
(red dot on left) chases a prey (red dot on right), while the 
latter tries to escape through chemotactic gradient sensing of 
the diffusing chemicals. The arrows indicate their respective 
direction of motion in absence of fluctuations. The contours 
around each microbe represent the equi-concentration lines of 
the secreted chemicals in a two-dimensional projected plane 
in this case, indicating the asymmetry of the distribution. 
The color code used here for the spatial distribution of the 
secreted chemorepellant (ci) and the chemoattractant (02), 
as they mingle in space, is shown in the right panel. 



diffusivity 300/im 2 /s. Microglial cells [26| moving at 
2/xm/min in response to an interleukin gradient of 

0. 003nM/^m secreted at 200 molecules/min and diffus- 
ing at 900/mi 2 /min, has ~ fOeo^g. 

A first look at the problem suggests that the ulti- 
mate fate of the prey, i.e. whether it will manage to es- 
cape, get captured or be steadily hunted forever, will 
depend not only on the emission rates and diffusivi- 
ties of the chemicals, the effective coupling strengths 
of the microbes to the chemical and their mobilities in 
the medium, but also on the initial separation between 
them: r = r 12 (t = 0) = |r 2 (0) — i"i(0)|. Having set the 
model, it is therefore instructive to examine the zero- 
noise (?7j(i) = 0) deterministic case first, in order to 
understand the combination of variables relevant in pre- 
dicting the outcome. Fig.l shows a simulation snapshot 
of the chemotactic chase process in the absence of fluctu- 
ations. Assuming a steady-state velocity Vi for microbe 

1, moving along the rc-axis, the concentration profile for 
the chemical secreted by it simplifies to 

d(r,t = 0) = (Xi/27TD ci r)exp(-Vi(x + r)/2D ci ) (7) 

In the presence of chemotactic coupling, one then ex- 
pects jiVi = (1 — 5ij)Ki [dej (r) /dx] T — Ti , , for each of them. 
In addition, demanding a condition for steady hunt- 
ing, whence Vi = Vj, maintaining a constant separation 
ri2 = A, leads to 



5(A*) = (1 + A*" 1 )exp(-A* 



(8) 



Here, A* = A/A , A = kiA 2 /(47t£) c iD C 27i) being a 
length scale, and 5 = (ki72A2-D c i)/(k27iAi-C> C 2) will be 
termed as a sensibility ratio in the predator-prey rela- 
tionship. Fig. 2(a) shows the resulting phase diagram. 
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FIG. 2: (color online). (a) Dynamical Phase Diagram of the 
chemotactic predator-prey system, constructed in the A*-S 
parameter space, showing the trapped (shaded) and escaped 
phases. The phase boundary (thick solid line) is obtained 
analytically, and matches with the simulation data (boxes). 
The horizontal thin dotted line (8 — 1) represents the upper- 
bound for the trapped-to-escaped dynamical phase transition 
(see text), (b) The dependence of the catching range (A*) on 
the initial separation (ro*), as obtained from simulations. 
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FIG. 3: (color online). Mean square displacement of the prey 
w.r.t. the predator as a function of time, with Ai = 1, 71 = 
10, 72 = 0.01, Dd = 10 3 , D c2 = 10 2 , ki = 10 4 for (a) 
/3" 1 = 0.01,0.1, 1.0 and k 2 = 15.5 x 10 4 ; and (b) zero noise 
case (p- 1 = 0) and k 2 = 15.5 x 10 4 ,5.5 x 10 4 , 1.5 x 10 4 . 
The power law behaviors t 2 ^ 3 and t are illustrated by the 
corresponding reference lines drawn. 



For a given A*, if the sensibility ratio is increased, there 
is a transition from the free to the trapped phase. To 
understand this we note that an increase in the chemo- 
attractant coupling or its emission or a decrease of the 
prey's mobility would prove advantageous to the predator 
in sensing the prey at a given distance. Also, a decrease 
in the chemoattractant diffusivity will enable the preda- 
tor to easily track its prey and ultimately trap it. Simi- 
larly, for a given S, increase in the separation distance will 
be advantageous to the prey in escaping. Further, since 
A > (predator follows prey), the validity of Eq.(8) re- 
quires 6 < 1. This means the phase boundary between 
the trapped and escaped state lies below the S(A*) = 1 
line in the parameter space. For 8 > 1, there is no es- 
cape. Whatever is the separation between the microbes, 
the predator will ultimately capture the prey in this case. 

In our simulations, the control parameter is the ini- 
tial distance ro between predator and prey. By tuning 
5 for a given ro, we actually found the point of tran- 
sition from trapped to escaped state. The border-line 
case of steady hunting was found to be unstable. Close 
to the border-line situation, the predator-prey distance 
ri2 remains constant for long times before fluctuations 
throw them into either the trapped or the escaped state. 
This constant distance, if identified with A, matches the 
phase boundary (Eq.(8)) perfectly, and A* depends on 
r o(= r o/Ao) in a roughly linear fashion (Fig. 2(b)). 

What are the dynamical features of the escaped and 
the trapped phases? Our simulations show that for es- 
cape (Fig. 3(a)), the mean squared displacement of the 
prey w.r.t. the predator gradually deviates from the 
initial value ro and grows subdiffusively with time as 
ri2 2 ~ t a , with an exponent a = 2/3, where the bar 
denotes an averaging over noise for a given ro. This be- 
havior finally crosses over to diffusion, ri2 2 ~ t, for long 
times. The crossover time t co decreases with increasing 
fluctuation strength Zero-noise simulations show 

the subdiffusive motion as the final long-time behavior, 
with the same exponent a (Fig.3(b)). This means that 



within this phase, the hunting process continues with a 
subdiffusive dynamics of the prey in the comoving frame 
of the predator; but finally due to effects of fluctuation, 
the prey diffuses away freely. It is therefore appropriate 
to look for a theoretical estimate of a within the noise-less 
ideal case. We note from Eq.(7) that at the advancing 
predator position, the chemo-attractant profile behind 
the prey is of the form C2(|x|) = A2/ (4:irD C 2\x\) in the 
steady state condition, since x < 0. Solving the resulting 
equation of motion in steady state, 71 ± = Ki|Vc2(x)|, 
gives x 2 ~ i 2 / 3 explaining the subdiffusive exponent. 
When the fluctuations ultimately overcome the chemo- 
tactic coupling, crossover to final diffusion results, re- 
quiring /3" 1 = KiC2(|x|) = Ki\ 2 /(4:TrD C 2\x\). At the 
crossover point from subdiffusive to the diffusive regime, 
x(t = tco) ~ tco 1 / 3 ; implying that the crossover time 
scales with the inverse effective temperature as t co ~ ft 3 ■ 

At the front of each microbe, the respective chemi- 
cal profile decays much faster ~ - exp (— Vix/D C i). The 
predator is thus always at an advantage of sensing the 
prey from much longer distances. Therefore, for very 
low S, a small increase in 6 greatly increases the catch- 
ing range for the predator. This accounts for the almost 
vanishing slope of the phase boundary for low 6 in the 
dynamical phase diagram (Fig. 2(a)). In this part of the 
phase diagram, for two predators with close sensibility 
ratio, the one with the slightly larger 5 will successfully 
trap preys which started at a much larger initial sepa- 
ration. For intermediate values of <5, the catching range 
A* increases at a slower rate with increase in S. This is 
because the initial separation is large enough that small 
changes in A* do not appreciably increase the chemoat- 
tractant gradient (~ 1/x 2 ) for the predator. One then 
needs to considerably alter the sensibility ratio for ob- 
taining a significant change in the chemotactic coupling. 
As (5 is further increased, A* increases to diverge at 5 = 1. 

For the trapping situation we find a broad skewed dis- 
tribution (P(tca P )) of capture time t ca p, for a fixed j3, tq 
and S (Fig.4(a)). The mean capture time t cap is, however, 
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FIG. 4: (color online), (a) Probability distribution of the 
capture time for the typical values of S and A* inside the 
trapped phase, for /3 = 100, 500, 2000. (b) Mean capture time 
(tcap) and variance (<r t ) of the capture time, as a function of 
/3. (c) Mean square predator-prey separation (ri2 2 ), close to 
trapping situations, as a function of the time interval |r| = 
\t — t cap \. The thick reference line indicates r 2//3 power- law 
behavior, (d) Divergence of the capture time with decrease in 
S, for a fixed catching range (A*), as the trapped-to-escape 
transition (S = S c , shown by the vertical line) is approached. 



independent of the effective fluctuation strength, while 
the variance at 



tcap) 2 ] 1 ^ 2 decreases with in- 



creasing (3 (Fig.4(b)). The trapping dynamics also show 
non-trivial power-law behavior, very close to capture: 
r 12 2 (r^0-) ~ |t| 2 / 3 , where r = t - t cap (Fig.4(c)). 
We note that initially, close to t = 0, when the predator- 
prey distance is ~ r$, fluctuations dominantly control 
the individual microbe's motion until a steady chemical 
concentration profile sets up in the process to favour a 



systematic dynamics. This time scale is dependent on 
the individual diffusivities of the chemical: higher the 
chemical diffusivities compared to the non-chemotactic 
diffusivities of the microbes, the faster will the predator- 
prey chemotactic systematics set in. For the trapping 
dynamics very close to capture, therefore, the predator 
is already responding to the steady chemoattractant gra- 
dient 7ii = Ki \Wc2(r)\ r — x ~ l/x 2 . Thus, integrating 
from x = x(t) to x — x(t — t cap ) = 0, we obtain the scal- 
ing form x 2 ~ (t cap — i) 2 / 3 , with t < t cap . For fixed initial 
separation, t cap diverges on approaching from above the 
corresponding critical value of the sensibility ratio (<5 C ) 
for escape (Fig. 4(d)). 

In conclusion, we studied the dynamics of simple gra- 
dient sensing chemotactic microorganisms in a predator- 
prey relationship. Although prokaryotes like most bac- 
teria chemotax chiefly through temporal comparison of 
chemical gradients, the direct spatial sensing which we 
studied here is prevalent among eukaryotes like amoeba, 
yeast cells, neutrophils, lymphocytes and glial cells. Here 
we have been able to delineate power law behaviors 
in the chasing process both from theory and numerical 
simulations, and crossover time scaling with fluctuation 
strength. A dynamical phase diagram has been obtained 
to identify conditions for escape and catching, with a 
border-line unstable steady-hunting situation. A broad 
class of microbes varying widely with respect to their mo- 
bility, secretion rates and diffusivities of ejected chemicals 
and strength of spatial gradient sensing of chemicals can 
be located in the phase diagram. Interestingly, a sim- 
ple sensibility ratio resulting from our model calculation 
allows a simple criterion to predict the outcome of such 
hunting processes: a trapped or an escaped situation de- 
pending on the initial predator-prey separation, and a no 
escape situation independent of their initial separation. 
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